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Abstract 

We analyze nuclear matrix elements (NME) of neutrinoless double beta decay calculated for the Cadmium isotopes. 
Energy density functional methods including beyond mean field effects such as symmetry restoration and shape mixing 
O^l are used. Strong shell effects are found associated to the underlying nuclear structure of the initial and final nuclei. Fur- 
thermore, we show that NME for two-neutrino double beta decay evaluated in the closure approximation, M^" , display 



O 



a constant proportionality with respect to the Gamow- Teller part of the neutrinoless NME, Mq^. This opens the pos- 
sibility of determining the Mq^ matrix elements from /3 T Gamow- Teller strength functions. Finally, the interconnected 
role of deformation, pairing, configuration mixing and shell effects in the NMEs is discussed. 
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1. Introduction 

Neutrinoless double beta decay (0vf3/3) is the most 
promising process to disentangle the Majorana nature of 
the neutrino, its effective mass and the mass hierarchy [f]. 
In this process, an even-even nucleus can not energetically 
decay into the odd-odd neighbor but it is allowed into the 
even-even nucleus with two protons more and two neutron 
less. Only few candidates which fulfill this requirement 
are found along the nuclear chart in the valley of the sta- 
bility. Contrary to the two-neutrino double beta decay 
mode (2vf3f3), where two electrons and two Dirac or Ma- 
jorana neutrinos are emitted in the final state, the 0^/3/3 
mode proceeds by exchanging a Majorana neutrino and 
only two electrons are emitted at a sharp energy equal to 
the Q-value. 2z//3/3 decay has been observed in several iso- 
topes with half-lives I0 19-21 years. However, except to 
the controversial claim of the Heidelberg-Moscow exper- 
iment [2], the neutrinoless mode has not been measured 
yet. Nowadays, several experiments devoted to detect this 
decay mode are presently active or in an advance develop- 
ment phase [3, 4]. If this process is finally detected, the 
precise determination of the neutrino effective mass de- 
pends on the value of the nuclear matrix element (NME) . 
The half-life of this decay, in the so-called light neutrino 
exchange mechanism, can be written as [1]: 
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(mpp) is the effective Majorana neutrino mass. Several 
nuclear structure methods have been used so far to cal- 
culate these NMEs for the most promising candidates, 
namely, large scale shell model (SM) [6, 7], quasiparti- 
cle random phase approximation (QRPA) [8-10], inter- 
acting boson model (IBM) [11], projected Hartree-Fock- 
Bogoliubov (PHFB) [12] and energy density functional 
methods (EDF) [13, 14]. In the most recent calculations, 
the spread in the value of NMEs between different models 
for a specific candidate is about a factor two. Therefore, it 
is necessary to study in more detail the NMEs in order to 
better constraint their values and understand possible re- 
lationships between the magnitude of the matrix element 
and the nuclear structure of parent and daughter states. 
In this Letter we use energy density functional methods 
including beyond mean field effects to analyze the NMEs 
for the decay of Cadmium to Tin isotopes. All of these 
decays except the one of 116 Cd are not physically possible 
because of the Q-value or because they are strongly hin- 
dered by single (3 decays. However, these virtual Qvj3f3 
decays provide useful information about the dependence 
of the NMEs with the underlying nuclear structure of the 
mother and daughter nuclei. Therefore, shell effects or the 
role of deformation and pairing can be studied in a system- 
atic and controllable manner. The paper is organized as 
follows. The theoretical formalism is described in Sec. 2. 
Then, the results obtained are reported in Sec. 3. Finally, 
a summary of the main conclusions is given in Sec. 4. 



where Gq u is a well-known kinematic phase space fac- 
tor [5], M 0v is the NME, m e is the electron mass and 
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2. Theoretical framework 

We now describe the main aspects of the generator 
coordinate method (GCM) used to compute the NMEs 
with energy density functional methods. In this frame- 
work, we determine the initial (i) and final (/) states as 
linear combinations of particle number and angular mo- 
mentum projected Hartree-Fock-Bogoliubov (HFB) mean- 
field wave functions [15-19]: 
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Then, the states given in Eq. 4 can be expressed as |-rtp = 
Ea 4/ G A z/f \ A l/f) and the HGW equations (Eq. 3) read 
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The eigenvectors G\. f are used to define the relative weight 
of each deformation j3 2 in the GCM wave function, the so- 
called collective wave functions [15]: 
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where / is the angular momentum, (3 2 are intrinsic ax- 
ial quadrupole deformations, gljf(P 2 ) are the coefficients 
found by solving the associated Hill- Wheeler- Griffin (HWG) 
equations [15]: 
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Here, E^J^ are the energy spectra of the initial/final nuclei 
and the projected wave functions are defined as: 
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with P N , P z , P 1 being the projection operators onto 
good number of neutrons, protons and angular momen- 
tum, respectively. Furthermore, for each deformation p 2 , 
the HFB-type states \<f)(/3 2 )) are found by minimizing the 
particle number projected energy with constraints in the 
mean value of the axial quadrupole moment operator Q 2 q, 
i.e., 6(E' NZ (\0(p 2 )))) = with [17]: 



e 'nz = W2)\HP N P z \m)) 



(5) 



The Lagrange multiplier Xp 2 ensures the condition 
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with r Q = 1.2 fm and A the mass number. We use the 
same interaction both in the determination of the intrinsic 
states and in the HWG diagonalization, in our case the 
Gogny D1S force [20]. To solve the HWG equations we 
transform first Eq. 2 in terms of an orthonormal basis: 
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where u A .. (fiz) and nj^. are the eigenvectors and eigen- 
values -greater than a small value e ~ 10~ 5 chosen to re- 
move linear dependence of the states- of the norm overlap 
matrix: 

E<*i//(^)i^//^)X//(^) = ( 8 ) 



In addition, any expectation value between GCM wave 
functions (energies, radii, transitions, etc.) can be ex- 
pressed in terms of these coefficients, as we will see in 
the case of 0isf3/3 NMEs. 

It is important to underline that, in this framework, 
particle number and rotational symmetries are conserved 
by projecting the intrinsic HFB wave functions and, in ad- 
dition, shape mixing is naturally included and self-consistcntly 
determined by solving the HWG equations. On the other 
hand, the main restrictions of the model come from the 
symmetries imposed to the HFB-type states \4>(f3 2 )). They 
are tensorial products of proton and neutron wave func- 
tions (no isospin mixing), and are also parity and axially 
symmetric. Work is in progress to check the effect of these 
approximations in the NMEs. Finally, a large configu- 
ration space including eleven major harmonic oscillator 
shells and sets of 56 states with different deformation rang- 
ing from /3 2 G [—0.5,0.6] are used in this work. 

Once the initial and final states are found, we proceed 
to evaluate the NMEs as the sum of Fermi (F) and Gamow- 
Teller (GT) terms [1] (tensor contribution is neglected in 
this work [7, 10]): 
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with gy = 1 and g& = 1.25 being the vector and axial cou- 
pling constants. First, we use the closure approximation 
to by-pass the calculation of the odd-odd intermediate nu- 
cleus. This approximation is expected to be a good one in 
the 0i>i3j3 case [1, 21]. Then, the NMEs can be determined 
as the expectation value of two-body operators between 
the initial and final states: 
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In these expressions, f_ is the isospin ladder operator that 
changes neutrons into protons and <x are the Pauli ma- 
trices acting on the spin part of the wave functions. The 



so-called neutrino potentials Vf/gt depend on the relative 
distance between two nucleons. These potentials take into 
account nucleon finite size corrections and higher order 
currents [8] (see Refs. [7, 8] for detailed expressions). Be- 
fore including short-range correlations within the UCOM 
framework [9, 22] the neutrino potentials can be expressed 
as an integral in the momentum transferred q of the Fermi 
or Gamow- Teller form factor hp/Qxil) weighted by the 
spherical Bessel function jo{qr) [7, 8]: 
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For the shake of simplicity we take the same value of the 
closure energy // = 10.22 MeV for the whole isotopic chain. 
This value is the one used for computing the NME of 116 Cd 
decay in Ref. [13]. 

We now calculate the NMEs (Eq. 12) with the EDF 
formalism described above. We take the ground states 
given by the lowest energy solution of the HWG equation 
with (I = 0, a = 1) and evaluate: 
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where M^ GT (/3 2 , /3 2 ) = (^ 2 )\M°/ /GT \^((3' 2 )). The 
last expression, properly normalized, gives the explicit de- 
pendence of the NMEs with the quadrupole deformation 
of the mother and daughter nuclei: 
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(17) 

The final value of the NME can be approximately inter- 
preted as the convolution of the collective wave functions 
of initial and final states (Eq. 10) with the intensity of 
the Fermi and GT NMEs as a function of the deformation 
given in Eq. 17. 

3. Results 

We now show the results of the 0v/3/3 NMEs in the 
cadmium isotopic chain 98_132 Cd. This set of nuclei covers 
the shells from iV = 50 to N = 82 magic numbers both 
in the mother (Cd) and daughter (Sn) isotopes. We first 
describe the ground states of these nuclei that are the ones 
involved in the Qv(3f3 decay. In Fig. 1 the experimental [23] 
and theoretical Qp-p-- values are plotted. We observe 
that experimental and theoretical data almost run paral- 
lel each other, increasing the value continuously with in- 
creasing the mass number from negative to positive values. 
Additionally, a quantitative agreement is found within the 
limit of accuracy of the Gogny D1S force, which is approx- 
imately 3.4 MeV [24]. Nevertheless, the theoretical result 
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Figure 1: (color online) Experimental [23] and theoretical Q^_^_ 
value for cadmium isotopes as a function of the mass number A. 



overestimates the Qp-p- -value about ~ 1.4 MeV. It is im- 
portant to point out that no tuning of the parameters of 
the interaction has been made. 

In Fig. 2 we represent the evolution from shell to shell 
of the ground state collective wave functions (Eq. 10) for 
Cd -Fig. 2(a)- and Sn -Fig. 2(b)- isotopes. First of all, we 
observe that the distributions are enclosed around /3 2 ~ 
[—0.2,0.2] for all nuclei, i.e., these are not very deformed 
systems. In addition, most of the collective wave functions 
present two maxima, one oblate and one prolate, with a 
minimum in the spherical point. This is not the case for 
the magic or their nearest neighbor nuclei ( 98 > 13 °- 132 Cd 
and 98 - 102 ! 13 °- 132 Sn) whose collective wave functions also 
peak in the spherical point. For the Cd chain (Fig. 2(a)) 
prolate maxima are higher than the oblate ones while a 
more symmetric distribution around j3 2 = is obtained 
for the Sn isotopes. In addition, mirror nuclei 98 Cd and 
98 Sn have almost identical collective wave functions as one 
could expect. 

A more quantitative analysis of the deformation can 
be done by representing in Fig. 3(a) the mean value f3 2 
for the collective wave functions given above. We see in 
Fig. 3(a) that tin isotopes are spherical or very little oblate 
deformed in the mid-shell (fi 2 ~ —0.025 for 114 Sn) as it is 
expected from its proton shell closure Z = 50. On the 
other hand, the mean deformation in cadmium isotopes 
is zero only in the shell closures V = 50, 82 while the 
rest of nuclei are slightly prolate deformed. Hence, we 
find two maxima at 104 . 106 Cd and 118 Cd with /3 2 - 0.12 
and f3 2 ~ 0.10 respectively and a local minimum at 112 Cd 
2 ~ 0.08). 

The effect of the deformation of initial and final states 
and the role played by pairing correlations on the NMEs 
for individual decays have been already discussed exten- 
sively within the EDF framework [13, 14], and, to a lesser 
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Figure 2: (color online) Ground state + collective wave functions 
for (a) cadmium and (b) tin isotopes calculated with the Gogny D1S 
EDF. 116 Cd ( 116 Sn) double beta emitter is represented in red. 



extent, within the shell model [25], QRPA [26] and IBM [11] 
approaches. On the one hand, a strong suppression of the 
NMEs is found for decays between states with different 
intrinsic deformations. Furthermore, shape mixing effects 
tend to the reduce the values of the matrix elements calcu- 
lated assuming only spherical symmetry [14, 27]. On the 
other hand, the NMEs are enhanced if a large amount of 
pairing correlations are present, or in the shell model lan- 
guage, when the wave functions are dominated by general- 
ize seniority zero components, in the mother and daughter 
nuclei [6, 13]. Let us study now if these conclusions can be 
extended to the NMEs calculated in the Cd isotopic chain 
where the number of neutrons changes smoothly within a 
major shell. 

In Fig. 3(b) we compare the values of the NMEs ob- 
tained from the full configuration mixing (Eq. 16) with 
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Figure 3: (color online) (a) Mean value of the quadrupole deforma- 
tion /?2 for the collective wave functions given in Fig. 2. (b) Nuclear 
matrix elements for initial and final nuclei considered as spherical 
systems (filled squares) and including the full shape mixing (bul- 
lets), (c) Total pairing energies for Cd and Sn isotopes including 
shape configuration mixing. 



those assuming spherical shapes M "^ — 0, (3' 2 = 0) 
(Eq. 17) for mother and daughter nuclei. We observe that, 
except for the transition between mirror nuclei A = 98, 
M 0l/ is small near the neutron shell closures. Further- 
more, the NMEs for the spherical shape are always larger 
than the ones which include configuration mixing although 
both curves show a similar structure. The reduction ob- 
tained by shape mixing is smaller near the shell closures 
where the relevant deformations explored by the collec- 
tive wave functions are precisely the spherical ones. From 
Fig 3(a) we can also see qualitatively that the larger is 
the difference between the mean quadrupole deformation 
of mother and daughter nuclei the larger is the reduction 
of the NME with respect to its spherical value. This shows 
that the correlation observed between the value of NME 
and differences in nuclear deformation is a general feature 
of neutrinoless double beta-decay and not a particular as- 
pect of those beta-decays studied so far [13]. 
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Figure 4: (color online) Nuclear matrix elements calculated with 
Fermi (squares), Gamow- Teller (bullets) and constant neutrino po- 
tentials Vf{t) = c = 2.0 (boxes) and Vgt(t) = d = 2.0 (circles) 

We now analyze the structure of the NMEs along the 
isotopic chain (Fig. 3(b)). First, we observe a large value 
of the NME for A — 98 transition between mirror nuclei. 
The wave functions of the initial and final states are prac- 
tically identical (see Fig. 2) and their overlap is therefore 
maximized [25]. Apart from this specific case, the rest of 
the curve shows some structure that can be related to neu- 
tron sub-shell closures in Cd isotopes. In Ref. [11] shell 
effects are studied within a generalized seniority scheme 
and it is observed that the lowest NME values correspond 
to the shell closures and a maximum is found at the mid- 
shell. The same behavior can be also seen in Fig. 3(b) 
where three maxima are found at A(N) — 102(54), 112(64) 
and 122(74) that can be related to the half-filling of 2d 5 / 2 , 
lff7/23si/22d 3 /2 and l/iu/2 sub-shells. The local minimum 
found at A = 118 would correspond to the clossing of the 
i57/23si/ 2 2d 3 /2 sub-shell. 

These shell effects are also related to the content of 
pairing correlations in initial and final nuclei as we can 
see in Fig 3(c). In this figure we represent the pairing 
energy (—E pp ) [15] evaluated with the GCM states given 
in Eq. 2. In closed shell nuclei, N = 50 and N = 82, 
pairing correlations are significantly smaller than in open- 
shell nuclei, obtaining the parabolic shape given in the 
figure. Nevertheless, we also find local minima in the Cd 
isotopic chain at A = 104 and A — 118 related to the 
sub-shell closures described above. Therefore, we obtain 
again a direct generalization of a result already found in 
individual decays, namely, the utter connection between 
pairing correlations in the initial and final nuclei and the 
NMEs. 

Finally we study the sensitivity of the NMEs to the 
form of the neutrino potentials Vp/qt{t) given in Eq. 15. 
In Fig. 4 we compare the Fermi and GT components of 
M 0u assuming a constant neutrino potential, i.e. indepen- 



dent of the relative distance of the decaying neutrons, with 
those obtained using the usual form (15). We take as 2.0 
the constant value of the neutrino potential as it repro- 
duces both the Fermi and Gamow- Teller NME computed 
using the full potential. The same pattern is found for 
the NMEs calculated with a constant spatial dependence 
and the ones given by the full neutrino potentials, i.e., 
similar enhancement of M 0l/ in mirror nuclei and shell ef- 
fects both for Fermi and Gamow- Teller parts. This shows 
that the overall behavior of the NME is quite insensitive 
to the details of the neutrino potential. Nevertheless, the 
NMEs computed with regular neutrino potentials are more 
symmetric with respect to the local minimum found at 
A = 118, while the ones computed with constant poten- 
tials are smaller after this minimum. This fact is related 
to the filling of the negative parity l/i-ii/2 sub-shell and 
could be associated to a increased contribution of multi- 
poles of higher order that the Gamow- Teller, i.e. the only 
one present in the constant potential case. 

This proportionality is important because the GT part 
calculated with a constant potential is directly related to 
the 2j//3/3 matrix element evaluated in the closure approx- 
imation [21]: 

M cT = <0;i Yl ■ ^ U) r ( - ] T U} |0+) (18) 

i<j 

Limits on this matrix element can be experimentally ob- 
tained from charge exchange reactions (see [28] and ref- 
erences therein) or assuming single-state dominance hy- 
pothesis [29]. Hence, it would help to constrain experi- 
mentally the value of M(V T . In this work we obtain for 
116 Cd an unquenched value of — 1.885 which is much 
larger than the one extracted from charge exchange re- 
actions [28] M^ ceI = 0.314. However, the GT operator 
would require a quenching factor to account for renormal- 
ization effects of the weak-axial current. On the other 
hand, the experimental data reported on [28] includes GT 
strength located at excitation energies below 3 MeV in 
the intermedium nucleus. One expects substantial con- 
tributions to the from states in the GT resonance 
at excitation energies around 10 MeV. Notice, that these 
transitions are suppressed for the regular M 2v matrix el- 
ement [21] due to the energy denominator but not for the 
M^i matrix element. Based on QRPA calculations ref. [21] 
shows that a converged value of M^ 1 requires the inclusion 
of states up to 15 MeV excitation energy in the intermedia 
nucleus. Further work is in progress to determine the ori- 
gin of this proportionality which has also been observed, 
within the EDF framework, in the NMEs as a function of 
deformation of individual decays [14], and in shell model 
calculations [25], but not found in QRPA calculations [21]. 

4. Summary 

In summary, we have studied the Qv(3f3 nuclear ma- 
trix elements in the cadmium isotopic chain with energy 
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density functional methods. It is shown that previous cor- 
relations found for individual decays between the magni- 
tude of the NME and pairing correlations, deformation, 
shell effects and shape of the neutrino potential can be 
extended to a whole isotopic chain. This shows that this 
correlations are a particular feature of the NME that is 
independent of the detailed structure of the nuclear wave 
functions. Large NMEs are obtained if the pairing corre- 
lations in the initial and final states are considerable and 
the difference in deformation is small. Furthermore, the 
NMEs are very sensitive to the sub-shell closures. In addi- 
tion, the NME for mirror nuclei is enhanced due to the 
large overlap between mother and daughter wave func- 
tions. Finally, NMEs calculated neglecting the neutrino 
potential, i.e. the two-neutrino double beta-decay matrix 
clement evaluated in the closure approximation, arc 1 found 
to exhibit a nearly constant proportionality to the realistic 
ones. This opens the possibility of constraining the M 0v 
matrix element from the knowledge of the Gamow- 
Tcller strength functions. 
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